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ANNOTATED  COMPUTER  OUPUT  FOR  SPLIT  PLOT  DESIGN:  GENSTAT  ANOVA 


W.T.  Federer,  Z.D.  Feng,  and  N.J.  Miles-McDermott 
Mathematical  Sciences  Institute,  Cornell  University,  Ithaca,  NY 

ABSTRACT 

In  order  to  provide  an  understanding  of  covariance  analyses  for 
split  plot  designs,  analyses  are  first  conducted  on  a  data  set 
without  the  covariate.  Then  an  analysis  with  the  covariate  is 
performed.  The  third  example  has  a  different  experiment  design  for 
the  whole  plots  and  the  covariate  is  constant  for  all  split  plots 
within  a  whole  plot.  Computer  packages  may  only  give  portions  of 
the  analysis  correctly.  Some  packages  require  two  procedural  calls 
to  obtain  a  portion  of  the  correct  results.  GENSTAT  requires  only 


one  procedural  call. 


INTRODUCTION 


This  is  part  of  a  continuing  project  that  produces  annotated 
computer  output  for  the  analysis  of  balanced  split  plot  experiments  ' 

with  covariates.  The  complete  project  will  involve  processing  three  (  ropv^ 
examples  on  SAS/GLM,  BMDP/2V,  SPSS-X/MANOVA ,  GENSTAT/ ANOVA ,  and 
SYSTAT/MGLH.  Only  univariate  results  are  considered.  We  show  here  Jj  . 
the  results  from  GENSTAT  ANOVA.  pj 

n 

For  Example  1,  the  data  are  artificial  and  were  constructed 
for  ease  of  computation;  the  experiment  design  for  the  whole  plots 

is  a  randomized  complete  block  and  the  split  plot  treatments  are  _ , 

1 1 

randomly  allocated  to  the  split  plot  experimental  units  within  each  - - 

whole  plot.  Example  2  is  the  same  as  Example  1  except  that  a 

fi-l , 

•.  •  <  .■  ,■  ,•  /  -■  „•  t, .-  v,  »•  .• w.  «■  f.  /.v  v,y,’<vV.> 


WWW 


w: 


covariate  varies  from  split  plot  to  split  plot.  The  data  for 
Example  3  come  from  an  experiment  wherein  the  whole  plot  treatments 
were  laid  out  in  a  completely  randomized  design  and  the  split  plot 
treatments  were  randomly  allotted  to  the  split  plot  experimental 
units  within  each  whole  plot.  The  value  of  the  covariate  varies 
from  whole  plot  to  whole  plot  but  is  constant  for  all  split  plots 
within  a  whole  plot  treatment. 

We  present  the  elementary  computational  steps.  Simple 
hypothetical  data  are  used  for  the  first  two  examples  so  that  it  is 
easy  to  provide  all  detailed  computations  to  illustrate  how  each 
number  is  obtained.  Some  readers  may  wish  to  skip  the  detailed 
computations.  The  third  example  comes  from  Winer  (1971).  The 
detailed  computations  are  given  in  his  book  (p.  803) . 

Data  SP-1 

Split  plot  data  with  whole  plots  arranged  in 
randomized  complete  block  design 
(hypothetical  data) 


W1 

Whole 

plot  treatment 

W2 

split 

plot 

treatment 

spl  it 

plot 

treatment 

Block 

S1 

S2 

S  3 

S4 

Total 

S1 

S2 

S3 

S4 

Total 

1 

3 

4 

7 

6 

20 

3 

2 

1 

14 

20 

2 

6 

10 

1 

11 

28 

8 

8 

2 

18 

36 

3 

6 

10 

4 

4 

24 

10 

8 

9 

13 

40 

Total 

15 

24 

12 

21 

72 

21 

18 

12 

45 

96 

1 


A  A  Vt  »>«  *«r.%‘iH  ^r»U‘ Ji  AVl 


a 


Blocks 

(8  observations) 


Total  and  Means 

W(whole  plots) 
(12  observations) 


S(split  plot) 

(6  observations) 


Total 


40 


Grand  Total 
Grand  Mean 


Mean 


5 


Total 


72 


Mean 


6 


Total 

36 


Mean 

6 


Model 


=  )J.  +  p  .  +  t  .  +  6.. 

*3  i  13 


+  €  .  . , 
i]k 


=  mean 


P3 
6  .  . 
13 

€ijk 


=  effect  of  block  j 
=  error  (a) 

=  error  (b) 


t.  =  effect  of  whole  plot  i 
a.  =  effect  of  split  plot  k 
(ar) =  effect  of  interaction  of 

whole  plot  i  and  split  plot  k 


Analysis  of  Variance 


Source 


B  (Blocks) 

W  (whole  plot  treatments) 

BxW  (error  (a)) 

S  (split  plot  treatments) 

SxW  (interaction  of  S  and  W) 

(**)  SxB:W  (error  (bn _ 

Total  (Corrected  for  mean) 
Mean 

Total  (Uncorrected  for  mean) 


R(p  m  ,  t  ,  a  ,  czr  ) 

R(t  u  ,p  ,a , ar ) 

R(<5  m  ,P  ,r  ,a,ar) 

R (a  M , P ,r  ,qt  ) 

R(ar lu ,a ,t , p) 

R  ( €  1  ji  ,  a  ,  t  ,  a  t  ,  p  ) _ 

R(p,T,6,a,ar,e |p) 
R(M) 

H  ill  ,p  J  T  ,  6  ,a  .ar  ,€  ) 


SS 


43 

24 

16 

156 

84 

_ 1 1 2_ 

440 

1176 

1616 


(*)Notation  follows  that  of  Searle ( 1971 ) ;  since  the  design  is  balanced, 
R(p  |p  ,  t  ,a  ,ar  )  =  R  ( p  |  jj  )  ,  etc.  The  simpler  notation  is  used  later. 

(**)  SxB:W  means  SxB  within  W. 

Calculations  of  sums  of  squares: 


N  =  2  •  3  •  4  =  24  ,  Y  =  7 

2  3  4 

R(U  ,P  ,t ,6,a ,ar ,€)=  III  Y^.  = 

i=l  j=l  k=l  ]K 


R(jj)  =  NY2  =  24*  (7)  2  =  1176 


2  2  2  2  2 
(3  +  6  +  6  +•  •  •  +  18  +  13  ) 


R(p ,t ,6 ,a ,ar ,e  |p) 


1616 


1176 


440 


R(p!m)  =  R(p,p)  -  R(p)  =  i40  +64+64  )  _  11?g  =  1224  _  1176  =  48 

2  2 

R(t|m)  =  R(u,T)  -  R(m)  =  i12^§-)  -  1176  =  1200  -  1176  =  24 

R(6|m,p,t)  =  R(6,4,p,t)  -  R(m,p)  -  R(t,m)  +  R(m) 

-  12_oW  +  2424202  +  362  +  402)  .  122„  .  120Q  +  ,176 
4 

=  1264  -  1224  -  1200  +  1176  =  16 

R(a|p)  =  R (a,u)  -  R(fu)  =  i^6  +42^  +  24  +66  )  _  11?g  =  1332  _  117g  =  156 

R  (aT  I M  / a  > T  )  =  R(ar,M,a,T)  -  R(p,a)  -  R(m,t)  +  R(p) 

=  122^l2421241824122t452)  .  1332  .  120Q  +  n76 

=  1440  -  1332  -  1200  +  1176  =  84 

R(t  |M,P(5,a,r,QT)  =  R(e ,u ,a ,p ,6 ,t ,ar) -  R(p,p,r,6)-  R(p,a,T,aT)  +  R(r,p 

=  1616  -  1264  -  1440  +  1200  =  112 


Data  SP-2 

Data  SP-2:  Data  SP-1  with  the  following  covariate  Z 
which  varies  with  split  plot 

Covariate  (Z) 


S1 

S  2 

whole  plot 

W1 

— 

— 

W2 

— 

S  3  S4 

Total 

S1 

S2 

S  3 

S4 

Total 

B1 

1 

2 

1  2 

6 

2 

0 

2 

4 

8 

B2 

2 

2 

0  4 

8 

4 

1 

3 

4 

12 

B3 

3 

5 

2  0 

10 

3 

2 

4 

7 

16 

Total 


2 


36 


Totals  and  Means 


blocks  W  (whole  plot)  S  (split  plot) 

(8  observations)  (12  observations)  (6  observations) 


Total 

Mean 

Total 

Mean 

Total 

Mean 

1 

14 

14/8  1 

24 

2.0 

1 

15 

2.5 

2 

20 

20/8  2 

36 

3 . 0 

2 

12 

2 . 0 

3 

26 

26/8 

3 

12 

2 . 0 

Grand 

Total 

60 

2 . 5 

4 

21 

3 . 5 

Model 

•  Y  .  = 

'  ijk 

u+p . +t . +6 . .+a 
3  l  13 

k+(«T>ik+/5 

1 v  ID  *  *  *  * 

)+/32 

'  13k  13- 

)+£ijk 

Pj 

°k 

=  effect 

of  jth  block 

t .  =  effect  of 
1 

ith  whole 

plot 

=  effect 

of  kth  split 

plot 

P1  =  whole 

plot  regression  slope 

^2 

=  split 

plot  regression  slope 

6^.  =  error 

a 

£ijk  ■  error  b 

Table  of  sura  of  squares  and  products 


Source 

df 

YY 

YZ 

ZZ 

B 

2 

48 

18 

9 

W 

1 

24 

12 

6 

Bxw  (error  a) 

2 

16 

4 

1 

s 

3 

156 

33 

9 

SxW 

3 

84 

33 

21 

SxB: W  (error  b) 

12 

112 

17 

20 

Mean 

1 

1176 

420 

150 

Total 

24 

1616 

537 

216 

YY  column  is  the  same  as  in  SP-1,  ZZ  column  is  computed  in  the  same 
fashion.  Thus,  only  computations  for  YZ  column  are  illustrated. 


T°talYZ 


2  3  4 

=  2  2  2  Y 

i=l  j=l  k=l 
-  3(1)+  6(2)+- 


i jk’ Zi  jk 

*•+  14(4)+  18(4)+  13(7) 


537 


Mean  „  =  NY  Z 

YZ  . 


168* 60 
24 


420 


W 


YZ 


3 

2 

=i zl 


2  4 

(  2  2 
i=l  k=l 


Yi  jk 


2 

)  (  ^ 
i  =  l 


4 

2 

k=l 


2*4 


2  3  4  3  4 

2(2  2  Y. ,  )(  2  2 

1=1  j=i  k=i  j=i  k=l 


3(4) 


420  -40[l4)+64  (20)^  64  (26) 
8 

=  438  -  420  =  18 


-  420  =  432  -  420  =  12 


420 


2  <2  Yijk><  *  Zijk> 

1  1  j — 1  K  1 - - - -  -  438  -  432  +  420 


=  454  -  438  -  432  +  420  =  4 

2  3  2  3 

4  *  Yijk><  2  .2  Zijk) 

g  .  j  —  1-^~  - - - - 1-1— -  -  420  =  453  -  420 

YZ  k=l  2(3) 

2  4  3  3 

2  2  (2  Y .  ..  )  (  2  Z .  ) 

i=l  k=l  i  =  i  ^  ^  k  i  =  l  ^k 

r-.,n  .  A  *  ^  A  J  J  a  rr  A  ^  i  A  r\ 


SxWYZ: 


-  453  -  432  +  420 


=  498  -  453  -  432  +  420  =  33 


SxB:Wyz: 


2  3  4 

2  2  2  Y...Z...  -  454  -  498  +  432 

i=l  j=l  k=l  13K  3 

=  537  -  454  -  498  +  432  =  17 


Analysis  of  Variance  and  Covariance 


B  (block) 

W  (whole  plot  treatment) 

=  R(P  U,r) 

=  R (t  ii,p,P1) 

Regression  (a) 

=  R(P1\p ,p ,t) 

B*W  (error  (a) ) 

=  R(5|m,P,t,P1) 

S  (split  plot  treatment) 

=  R(a  \p  ,  P  ,t  ,ar  ,  P  2) 

SxW  (interaction  of  S  and  W) 

=  R(ctT  \p  ,  p  ,t  ,a  ,P  2) 

Regression  (b) 

=  R(P2  Im  ,  P  ,T  ,a  ,clt) 

SxB:  W  (error  (b) ) 

=  R(fc |m ,p ,a ,t ,ar ,P  ) 

3 .4286 
16.0 
0.0 

84.243 
37 . 474 
14.450 
97 . 550 


*1  =  B 


The  SS's  adjusted  by  regression  on  Z  are  illustrated  below: 


R(P |m)  =  48,  remains  same  since  it  is  not  of  interest  to  adjust  for  Z 
on  the  blocks. 


R(t,6|u,p,/31)  =  (Wyy+  BxWyy)  - 


(Wyz+  B*Wyz) 


wzz+  BxWzz 


256 


=  (24  +  16)  -  -^1-2  +  4-^  =  40  -  =  3.4286 

6  +  1  7 


R(5|ja,p,T,P1)  =  BxW 


YY 


— .  16  .  ii  .  o 


BxW 


ZZ 


R(Tl4,P/P1)  =  R(T/6  I  M  ,  P  ,  /3  1 )  -  R(5|fi,p,T,/31) 


256 


=  40  -  -  0  =  3.4286 


.  (BXWYZ}  4 2 

R(P1|p,T,p)  =  - - -  =  —  =  16 


BxW 


ZZ 


R(a,e \u,p ,t ,ar ,p0)  =  (Syy  +  SxB:Wvv)  - 


yy> 


(SYZ  +  SXB:WYZ) 

szz  +  SxB:Wzz 


=  (156  +  112)  -  -£3-3-  17 ^ 

9  +  20 


=  268  -  86.207  =  181.793 


(SxW  +  S  x  B : W  ) 2 

R(cit  ,6\n,p,a,T,p  )  =  (SxW  +  Sx  B :  W  ) - — - — - 

SxWzz  +  SxB:Wzz 


=  84  +  112  -  ^-33  +  17^  =  196  -  60.976  =  135.024 

21  +  20 

Note:  R  (a  ,  e.  |p  ,  f?2 )  and  R(ar  ,e \p  ,a  ,t  ,p  are  intermediate  steps  for 

later  use. 


=  112 


=  112 


14.45 


97 . 55 


R(e\u  ,p  ,a  ,t  ,ar  ,P2)=  SxBiW^- 


(SxB:Wyz) 


SxB:  W 


ZZ 


17 

20 


R(a  |n  ,  p  ,t  ,ar  ,P2)  =  R  (a  ,  e  |p  ,  p  ,  t  ,ot  , /3  )  -  SS  error  b  =  181.793 

=  84.243 


-  97.55 


R(aT  Im  ,  P  ,a  ,T  ,P2)  =  R(ar  ,t\v  ,P  ,a  ,T  ,P2)  -  R  (t  |  ji  ,  p  ,  a  ,  t  ,  qt  ,  (5  2  ) 


=  135.024  -  97.55  =  37.474 
Data  SP-3 

Split  plot  data  with  plots  arranged  in  a  completely 
randomized  design  and  a  covariate  Z  that  is  constant 
within  the  whole  plot.  (Winer,  1971,  p.  803) 


whole  plot 

Subject 

Split 

B1 

plots 

B2 

j- 

! 

Cs3 

Total 

Y 

Y 

Y 

A1 

1 

10 

8 

3 

18 

2 

15 

12 

5 

27 

3 

20 

14 

8 

34 

4 

12 

6 

2 

18 

A_ 

2 

5 

15 

10 

1 

25 

6 

25 

20 

8 

45 

7 

20 

15 

10 

35 

8 

15 

10 

2 

25 

Total 

132 

95 

39 

227 

Mean 

16 . 5 

11.9 

00 

CO 

Model :  y...  =  u  +  t  .  +  6  .  . 

ilk  l  ii 

t ^  =  A  effect  (whole  plot) 
a,  =  B  effect  (split  plot) 


+  a.  +  (to).,  +  p .  ( Z  .  .  -  Z  )  +  t  .  ■ , 

k  v  ' lk  1 '  l  j  • • '  ilk 


<5 i j  =  error  (a; 


e  .  . ,  =  error  (b) 

ilk 


p.  =  whole  plot  regression  slope 


Analysis  of  variance  and  covariance 


Source _ df _ SS 


A  (whole  plot) 

=  R(t  |p  , /3 x ) 

1 

44 .492 

Regression 

=  RCPjJp/T) 

1 

166.577 

Error  (a) 

=  R(6  |m  ,t  ,/31) 

5 

61 . 298 

B  (split  plot) 

=  R(a |p ,t ,qt ) 

1 

85. 563 

AxB  (interaction) 

=  R(tq |p ,t ,a) 

1 

0 . 563 

Error  (b) 

=  R(e  lu  ,t  ,a  ,ra) 

6 

6 . 375 

Table 

of  SS  and 

products 

Symbol 

Y2 

ZY  Z2 

W 

68.06 

12.38  2.25 

E  (a) 

227.88 

163.00  159.50 

S 

85.563 

0  0 

WS 

0.563 

0  0 

-_E(b) _ 

6. 375 

0  0 

A 


163 . 00 
159.50 


1.02 


Since  the  computations  are  illustrated  in  Winer  (1971,  p.  803-5)  v 
have  omitted  them  here. 
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<5  i  c- 


Control  Language  is  typed  in  upper  case  and  cor.r.ents 


'REFE'  SPLIT 
'UNITS'  $  24 


FACTOR' 


file  name 

number  of  observations 
BLOCKS  $3 

\  define  levels  of  each  factor 


PLOTS  $2 
SUBPLOTS  $4 


1  READ/FORM=P, 

'RUN  ' 

1 

1 

1 

1 

3 

1 

1 

2 

2 

4 

1 

1 

3 

1 

7 

1 

1 

4 

2 

6 

1 

2 

1 

2 

3 

1 

2 

2 

0 

2 

1 

2 

3 

2 

1 

1 

2 

4 

4 

14 

2 

1 

1 

2 

6 

2 

1 

2 

2 

10 

2 

1 

3 

0 

1 

2 

1 

4 

4 

11 

2 

2 

1 

4 

8 

2 

2 

2 

1 

8 

2 

2 

3 

3 

2 

2 

2 

4 

4 

18 

3 

1 

1 

3 

6 

3 

1 

2 

5 

10 

3 

1 

3 

2 

4 

3 

1 

4 

0 

4 

3 

2 

1 

3 

10 

3 

2 

2 

2 

8 

3 

2 

3 

4 

9 

3 

2 

4 

7 

13 

BLOCKS,  PLOTS,  SUBPLOTS,  X,  Y 


I  nput  v ii  r  i  a  t 


’ EOD 1  =>  tell  GENSTAT  that  data  flow  ends 

'BLOCKS'  BLOCKS/PLOTS/SUBPLOTS  =r>  define  error  terms  (strata)  of  a 

factorial  model 

'TREATMENTS'  PLOTS *SUBPLOTS  treatment  terms  of  a  factorial  model 

'COVARIATES'  X  covariate 

'ANOVA'  Y  =>  invokes  analysis  of  variance  on  Y  variable 
'  RUN ' 

'STOP' 


SP-3:  Control  Language 


OK,  SLIST  SPLIT3 . OGEN 

•REFE'  SPLIT3 
'UNITS'  $  16 
'FACTOR'  PLOT  $2 

:  SUBPLOT  $2 
:  SUBJECT  $8 

' READ/FORM=P,  PRIN=D'  SUBJECT,  PLOT,  SUBPLOT,  Y,  Z 
'RUN' 

1  1  1  10  3 
112  8  3 

2  1  1  15  5 

2  1  2  12  5 

3  1  1  20  8 

3  1  2  14  8 

4  1  1  12  2 

4  12  6  2 

5  2  1  15  1 

5  2  2  10  1 

6  2  1  25  8 
6  2  2  20  8 

7  2  1  20  10 

7  2  2  15  10 

8  2  1  15  2 
8  2  2  10  2 
'  EOD ' 

'BLOCKS'  PLOT/SUBJECT/SUBPLOT 
'TREATMENT'  PLOT*SUBPLOT 
' COVARIATES '  Z 
' ANOVA '  Y 
'RUN' 

'STOP' 


SP-1:  Split  plots  with  whole  plots  arranged  in  RCB  design 
SP-2:  Split  plots  with  whole  plots  arranged  in  RCB  with  a 
covariate  with  split  plot 

Analysis  of  variance  table  on  covariate  Z 


*****  analysis  of  variance  ***** 


VARIATE :  Z 


SOURCE  OF  VARIATION 

DF 

SS 

SS% 

MS 

VR 

BLOCKS  STRATUM 

2 

9. 000 

13 . 64 

4 . 500 

BLOCKS. PLOTS  STRATUM 

PLOTS 

1 

6 . 000 

9 . 09 

6.000 

12 . 000 

RESIDUAL 

2 

1 . 000 

1 . 52 

0. 500 

TOTAL 

3 

7 . 000 

10.61 

2.333 

BLOCKS . PLOTS . SUBPLOTS 

STRATUM 

SUBPLOTS 

3 

9.000 

13 . 64 

3 . 000 

1 .800 

PLOTS . SUBPLOTS 

3 

21.000 

31.82 

7 . 000 

4 .200 

RESIDUAL 

12 

20.000 

30.30 

1.667 

TOTAL 

18 

50.000 

75.76 

2.778 

GRAND  TOTAL 

23 

66.000 

100.00 

GRAND  MEAN  2.50 

TOTAL  NUMBER  OF  OBSERVATIONS  24 


ANOVA  for  Y  variable  without  covariate  Z 


*****  ANALYSIS  OF  VARIANCE  ***** 

VARIATE:  Y 

SOURCE  OF  VARIATION  DF  SS  SS% 

MS  VR 

(F-statistics) 

BLOCKS  STRATUM  2  R  (p  |n  ,  t  ,a  ,clt  )  48.000  10.91 

24.000 

BLOCKS. PLOTS  STRATUM 


PLOTS 

1 

R(t \ii,p  ,a,ar)  24.000 

24 . 000 

5.45 

3 . 000 

RESIDUAL 

2 

R  (5  |n  ,  p  ,  t  ,  a  ,ar )  16.000 

8.000 

3 . 64 

TOTAL 

3 

40.000 

13.333 

9.09 

BLOCKS . PLOTS . SUBPLOTS 

STRATUM 

SUBPLOTS 

3 

R(a  |ji,p,T,aT,6)  156.000 

52.000 

35.45 

5.571 

PLOTS . SUBPLOTS 

3 

R(ar  ffx  ,a ,  T  ,  p  ,  S )  84.000 

28.000 

19.09 

3 . 000 

RESIDUAL 

12 

R(€  ,  a ,  t  ,ar  ,  p  ,  6)112. 000 

9.333 

25.45 

TOTAL 

18 

352 . 000 
19.556 

80.00 

GRAND  TOTAL 

23 

440.000 

100.00 

GRAND  MEAN  7.00 

TOTAL  NUMBER  OF  OBSERVATIONS  24 


*****  COVARIANCE  REGRESSIONS  ***** 


COVARIATE  (Z) 

COEFFICIENT 

SE 

BLOCKS  STRATUM 

Z  (PQ  is  not  of  interest) 

P0  "  2-0 

1 .  15 

BLOCKS. PLOTS  STRATUM 

BYZ  18 

Bzz  '  9 

Z 

A 

Pi  =  4 

0 

BLOCKS. PLOTS. SUBPLOTS  STRATUM 

Z 

P2  =  0.85 

0.666 

*****  TABLES  OF  MEANS  ***** 
(ADJUSTED  FOR  COVARIATE) 

VARIATE:  Y 

GRAND  MEAN  7.00 

PLOTS  1  2 

8.00  6.00  = 

Y2-.-  Pl(Z2..-  Z- 

..)  =  8-4. 0(3 

—2 . 5 ) =6 

SUBPLOTS  1  2 

3  4 

6.00  7.43 

4.43  10.15= 

Y..4-  ^2(Z-. 
11—0 .85 (— ^  - 

D 

11— . 85  =  10. 

.  -  Z 

A  •  •  • 

60 

24 ' 

15 

SUBPLOTS  1  2 

PLOTS 

1  7.00  9.15 

3  4 

6.85  9.00 

2  5.00  5.70 

2.00  11.30= 

Y2-4^1^2..- 

Z  ) 

•  •  •  ' 

*2-.> 

=  15-4(3  -  2.5)  - 

•  85 (— |  -3)  = 

11.30 

15 

f-V-  /‘V-V-jv' 

V„V  V  \r.v 

*****  STANDARD  ERRORS  OF  DIFFERENCES  OF  MEANS  ***** 


TABLE 

PLOTS 

SUBPLOTS 

PLOTS 

SUBPLOTS 

REP 

12 

6 

3 

SED 

0.000 

1 . 844 

2.718 

*****  STRATUM  STANDARD  ERRORS  AND  COEFFICIENTS  OF  VARIATION  ***** 


STRATUM 

DF 

SE 

CV% 

BLOCKS 

1 

1.225 

17 . 5 

BLOCKS. PLOTS 

1 

0.000 

0.0 

BLOCKS . PLOTS . SUBPLOTS 

11 

2.978 

42.5 

2.978  = 


unadjusted  standard  error 
■s/COV  EF  of  Residual 
V9.333 
Vl.052 


NOTE: 


S.E  (Y.  -  Y. ,  )  =  ( 

1*  *  1  *  *  l 


2E  W  /(W-l) 

a  (1  +  zz 


r  (s) 


BxW 


zz 


1  / : 


S.E  (Y.  Y  .  .) 


■{ 


0.0 

2E 


r(w) 


,  S  /(s-1) 

b  (1  +  zz 


Sx  B:  W 


z  z 


V: 


=  /  2(97.55)  +  g/(47lj  V/2  _ 

\  3 (2)  (11)  (1  20  } /  -  1.844 


S.E  (Y .  .  • -  Y  .  .  .  ) 
i  •  k  l  *  k 


■{ 


2E 


„  ,  (Szz+  SxWzz)/Ws-l)  v 

(1  +  - cTF.'u -  )  j 


SxB:  W 


zz 


=  /  1197^55)  19*21)7(2)  (3) 

\  3(11)  U  20  ;/ 


2.718 


where 


V4  = 
r  = 

Eb  = 


no.  of  whole  plot 
no.  of  blocks  =  3 


=  2 


S  =  no.  of  split  plot 


E  =  error  a  =  0 

a 


error  b 


97.55/11 


SP-3 


Split  plots  with  whole  plot  arranged  in  CRD  with  a  covariate 
constant  within  whole  plot 


*****  ANALYSIS  OF  VARIANCE  *****  ANOVA  on  Covariate  Z 


VARIATE:  Z 


SOURCE  OF  VARIATION 

DF 

SS 

SS% 

MS 

PLOT. SUBJECT  STRATUM 

PLOT 

1 

2 . 250E 

0 

1 .  39 

2 . 250E 

0 

RESIDUAL 

6 

1 . 595E 

2 

98 . 61 

2 . 658E 

1 

TOTAL 

7 

1 . 618E 

2 

100.00 

2 . 3  1 1 E 

1 

PLOT . SUBJECT . SUBPLOT 

STRATUM 

SUBPLOT 

1 

0. 000E 

0 

0 . 00 

0 . 000E 

0 

PLOT. SUBPLOT 

1 

0. 000E 

0 

0. 00 

0 . 000E 

0 

RESIDUAL 

6 

0.000E 

0 

0.00 

0 . 000E 

0 

TOTAL 

8 

0. 000E 

0 

0 . 00 

0. 000E 

0 

GRAND  TOTAL 

15 

1 . 618E 

2 

100.00 

GRAND  MEAN 

4  . 

88 

TOTAL  NUMBER  OF  OBSERVATIONS  16 

Notice  that  SS  in  plot . subject . subplot  stratum  are  all  zeroes 
since  Z  is  constant  within  whole  plot. 

ANOVA  on  Y  without  covariate  Z 

*****  ANALYSIS  OF  VARIANCE  ***** 

VARIATE:  Y 


SOURCE  OF  VARIATION 

DF 

SS 

SSI 

MS 

VR 

PLOT. SUBJECT  STRATUM 

PLOT 

1 

68.063 

17.52 

68 . 063 

1.792 

RESIDUAL 

6 

227.875 

58 . 66 

37 . 979 

TOTAL 

7 

295.938 

76.19 

42 .277 

PLOT . SUBJECT . SUBPLOT 

STRATUM 

SUBPLOT 

1 

85.563 

22.03 

85. 563 

80. 529 

PLOT. SUBPLOT 

1 

0. 563 

0. 14 

0. 563 

0.529 

RESIDUAL 

6 

6.375 

1 . 64 

1 . 063 

TOTAL 

8 

92 . 500 

23.81 

il . 563 

GRAND  TOTAL 

15 

388.438 

100. 00 

GRAND  MEAN 

14  . 

19 

TOTAL  NUMBER  OF  OBSERVATIONS 

16 

17 


Covariance  Analysis 


*****  tables  of  means  ***** 

(ADJUSTED  FOR  COVARIATE) 
VARIATE:  Y 


GRAND  MEAN  14 . 19 
PLOT  1  2 

12.51  15.87  =  V2#>-  P1(Z2.~  Z..)  =  16.25-1.02(5.25-4.88) 

=  15.87 

SUBPLOT  1  2 

16.50  11.88  same  as  unadjusted  mean 


SUBPLOT  1  2 

PLOT 

1  14.63  10.38 

2  18.37  13.37  =  *  2. 2~P  Z.  .  )  =  13.75-1.02  (5.25-4.88) 

=  13.37 

*****  STANDARD  ERRORS  OF  DIFFERENCES  OF  MEANS  ***** 

TABLE  PLOT  SUBPLOT  PLOT 

SUBPLOT 


REP  884 

SED  1.763  0.515  1.829 

EXCEPT  WHEN  COMPARING  MEANS  WITH  SAME  LEVEL ( S )  OF: 
PLOT  0.731 


*****  STRATUM  STANDARD  ERRORS  AND  COEFFICIENTS  OF  VARIATION  ***** 


STRATUM 

DF 

SE 

CV% 

PLOT. SUBJECT 

5 

2 .476 

17 . 5 

PLOT . SUBJECT . SUBPLOT 

6 

1.031  = 

Vl. 063 

Vl. 000 

7 . 3 

residual  MS  of  unadjusted  1 1/2 
.COV  EF  of  adjusted  residualj 


